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ABSTRACT 

iLutz & Kelkerl showed that parallax measurements are systematically overestimated because 
they do not properly account for the larger volume of space that is sampled at smaller par- 
allax values. We apply their analysis to neutron stars, incorporating the bias introduced by 
the intrinsic radio luminosity function and a realistic Galactic population model for neutron 
stars. We estimate the bias for all published neutron star parallax measurements and find that 
measurements with less than ~ 95% certainty, are likely to be significantly biased. Through 
inspection of historic parallax measurements, we confirm the described effects in optical and 
radio measurements, as well as in distance estimates based on interstellar dispersion mea- 
sures. The potential impact on future tests of relativistic gravity through pulsar timing and on 
X-ray-based estimates of neutron star radii is briefly discussed. 
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1 INTRODUCTION 



The past decade has seen an exponential increase in the number 
of parallax measurements to radio pulsars. Whereas only 14 such 
measurements were known by the end of the year 2000, currently 
52 radio pulsars have their distance determined either through 
VLBI or pulsar timing measurements (see Tables[TJand[2}. 

The importance of radio pulsar parallaxes arises from the 
wide variety of highly accurate investigations to which pulsars lend 
themselves. For example, the combination of an accurate paral- 
lax with the pulsar dispersion measure (DM: the integrated elec- 
tron density between the pulsar and Earth) provides an average 
electron density measurement along the line of sight, which can 
be used to construct Galactic electron density models (see, e.g., 
ICordes & Lazioll2002l) . Also, the highly polarised nature of pulsars 
allows measurement of the Faraday rotation which - in combina- 
tion with a distance - can be used to map out the co mponent of 
the G alactic magnetic field parallel to the line of sight Man et al.l 
2006). Finally, in pulsar timing, distances are essential in correct- 
ing spindown and orbital peri od derivatives for the Shklovskii ef- 
fect caused by proper motion lShkl ovskii|[T970h . Consequentially, 
certain pulsar timing tests of general relativity ar e dependent on ac- 
curate distance meas urements, as descr ibed by iDamour & Tavlotl 
dl99 lh and used in iNice et alj J2005t) and iDeller et alj J2008h . 
amongst others. 

The large increase in the number of radio pulsar parallaxes and 
their unique applicatio ns warrant an investiga tion into potential bi- 
ases. As pointed out bv lLutz & Kelken d 19730 . in an homogeneous 
field of stars the number of stars per unit of distance increases as 
D 2 , where D is distance. This makes it statistically more likely 
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to find an object at larger distance s where more v olume is sam- 
pled and, hence, more sources lie. ILutz & Kelkeil dl973l) showed 
analytically that this statistical underestimate of the stellar distance 
depends on the precision of the parallax measurement. Specifically, 
they derived the following proportionality: 



— 



exp 



(zu 



2a 2 



(1) 



where zuo is the measured parallax, a is the standard deviation of 
the measurement and p(zu\zuo) is the probabil ity distribution of the 
actual parallax, zu, given the measurement. iBinnev & Merrifiekj 
(l998t) expanded this analysis by using the intrinsic luminosity 
function for the star's spectral class as a further source of prior in- 
formation to use in the estimate of a more accurate parallax value. 

Two points of confusion have pervaded the literature on the 
topic of Lutz-Kelker bias. First, the effect we describe as Lutz- 
Kelker bias is commonl y referred to as M almquist bias in extra- 
galactic astronomy. As Gonzal ez & Faberl i 19971) point out, the 
Malmquist bias originally referred to a positive bias in luminosity 
of magnitude-limited samples and has onl y recently acqu i red th e 
meaning of th e geometric bias described bv lLutz & Kelken (1973). 
ISmithl J2003h stresses a second point of confusion and difference 
between the original Malmquist bi as and Lutz-Kelker bia s, namely 
the fact that the bias described by ILutz & Kelkerl i 19731) and dis- 
cussed in the present paper, refer to individual parallax measure- 
ments rather than to the overall average of a s ample of objects. 

In this paper we revisit the analysis by IBinnev & Merri field 

dl998h with a particular focus on radio pulsars. To account for the 
analytically complex but realistic pulsar luminosity function and 
Galactic pulsar distribution, we describe a Monte-Carlo approach 
to correct previously published parallax measurements for the Lutz- 
Kelker bias. Our basic analytic derivation and a description of the 
simulations are given in Sj2] The resulting corrections to published 
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measurements are discussed in Sj3] Our conclusions are summarised 
in g] 



2 THEORY AND SIMULATIONS 

As an illustration of the bias caused by volumetric and intrinsic 
luminosity arguments, it is helpful to c onsider an analytic de riva- 
tion along the lines of those presented bv lLutz & Kelkerf < fl973h and 
iBinnev & Merrifieldl dl998i) . Such an analysis is presented in £12 . 1 1 
below. For any practical purposes, however, a more realistic ap- 
proach that incorporates a Galactic distribution model for neutron 
stars as well as a more realistic intrinsic luminosity function for 
radio pulsars, is required. Because such realistic models are analyt- 
ically complex, we will implement them only as part of a Monte- 
Carlo simulation, which is described in §2.2\ 



2.1 Basic Theory 

For a pulsar at distance D — zu^ 1 , Bayes' theorem gives the prob- 
ability that the true parallax is raj, given the measurement zuq, as 



p(zu\zu ) = 



p(zu Q \zu)p(zu) 
p(ra?o) 



(2) 



where p(zuo) is the prior probability of the measurement which is 
taken to be constant without loss of generality. Assuming the mea- 
surement is dominated by Gaussian noise with standard deviation 
a, we know furthermore that 



2-Tva 



exp 



zu — Wq 



(3) 



The final term in Eq.[2]is the prior probability distribution of the 
parallax, p(zu), which is composed of two terms: one defined by the 
sampled volume and another defined by the intrinsic pulsar lumi- 
nosity function. Specifically, assuming a homogeneous distribution 
of pulsars in space, we have 

p(D) oc P D 2 , 

where p is the spatial density of pulsars. Translating this probability 
to a probability in m, we get 



dD 



dzu 



p(D) oc zu 



(4) 



Note this is effectively the prior lLutz&Kelkenjl973h determined. 
We will henceforth refer to this prior as the "volumetric" prior dis- 
tribution. 

The second contribution to the parallax prior is determined by 
the intrinsic radio pulsar luminosity function. Assuming, for sim- 
plicity, the power-law luminosity function 

p{L)^L p , 

where L is the intrinsic pulsar luminosity and j3 is the slope of the 
luminosity function, we obtain 



-2/3-3 



(5) 



which we will refer to as the "luminosity" prior distribution. Com- 
bination of Eqs.[4]and[5]results in the prior probability distribution 
of the pulsar parallax 



(zu) — p-o(zu)p-L(zu) oc w 



-2/3-7 



Inserting this along with Eq.[3]into Eq. [2] we obtain the resulting 
probability distribution of the true parallax: 

p(zu\zuo) oc zu~' 2l3 ~ 7 e 5 ( ° ) . (6) 

The (local) maximum of this probability distribution can easily be 
determined by requiring 

dp(zu\zu ) 



dzu 



0, 



which gives: 

(-2/3 -^^-VK^) 2 

o 

Rearranging this expression results in the analytic statement of the 
bias: 



zu 1 
~a ~ 2 



^W^-4(2/3 + 7) 
a V a 1 



(7) 



This result clearly demonstrates that the bias is a function of the 
significance of the parallax measurement and not of the parallax 
itself. 

Equation [7] also shows the effect of different intrinsic lumi- 
nosity functions: for steep power laws (J3 < —3.5) the luminos- 
ity function dominates the prior and the measured parallax value 
is likely to be underestimated. For shallow luminosity functions 
(/? > —3.5) the reverse is true and the parallax is expected to be 
overestimated. Note that in this case the prior in Equation [6] goes 
to infinity as zu goes to 0. A luminosity function with spectral in- 
dex ft = —3.5 results in a luminosity prior that exactly cancels out 
the volumetric prior, so that no bias is observed (as can easily be 
seen in Equ ation [6]). For the observed pulsar population, /3 ~ —1.7 
dLorimer et"afl l2006). which would imply measured parallaxes to 
be larger than their actual values. 

2.2 Realistic Monte-Carlo Simulation 

There are two clear improvements to be made to the analysis 
presen ted in the previous section. First, iFaucher-Giguere & Kaspi 
hood) have shown that the intrinsic luminosity function of radio 
pulsars is not a simple power-law but more likely log-normal in 
form. Second, pulsars are not distributed homogeneously in space, 
but are mostly confined to the Galactic disk. The analytic compli- 
cations implied by these improvements require that we proceed to 
evaluate the bias through Monte-Carlo simulations. Practically, we 
simulate the Gaussian measurement uncertainty and the two prior 
distributions independently, creating three normalised probability 
distribution^ Multiplication of these distributions provides the fi- 
nal, bias-corrected, distribution, p(zu\zuo) (see for example Figure 
[]}. We will report the peak of this combined probability distribution 
as the bias-corrected parallax value, zucon- 

In simulating the volumetric prior, we model a Galaxy of pul- 
sars wit h a radius of 15 kpc and with the radial density profile de- 
rived bv lLorimer et alj (2006) (their Equation 10, Model C fit): 



exp 



-C 



R — R, 



■© 



R 







1 Since both the scale height and intrinsic luminosity may be argued to 
depend on pulsar age, a more advanced, joint, simulation might be called 
for. Given current uncertainties in pulsar luminosities and the luminosity- 
age relation, however, we exclude this extension from the present analysis. 
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with constants A = 41 kpc -2 , B = 1.9, C = 5 and R = 8.5 kpc 
(i.e. the distance between the Sun and the Galactic centre). The 
density distribution above the Galactic plane is modelled by an ex- 
ponential distribution p(z) oc exp(—\z\/E) where E is the scale 
height, taken t o be 330 pc for comm on pulsars (in accordance with 
the findings oflLorimer et al.l2006|) and 500 p c for millisecond pul- 
sars dLorimeJl995l ; ICordes & Chernofil 19971) . To optimise the vol- 
umetric prior, we use the simulated Galaxy of pulsars to determine 
a volumetric prior probability density function for sky sectors of 
10° in Galactic latitude and 15° in Galactic longitude and for the 
polar regions, defined as having a Galactic latitude in excess of 85° . 

The prior information of the intrinsic radio pulsar lumi- 
nosity function is simulated through random realisations of pul- 
sar luminosities from the log- normal distribution derived by 
lFauche r-Giguere ~& Kaspil J2006I) . with mean intrinsic luminos- 
ity 0.07 mJy kpc 2 (i.e. (logi) = —1.1) and standard deviation 
o"io g L = 0.9 in base- 10 logarithms. This luminosity distribution 
is subsequently converted into a distribution of parallaxes using the 
known pulsar flux, S, via: Svj~ 2 = L. The flux values we used 
are given in column 3 of TableQ] Notice that, unlike the power-law 
luminosity function described in Section |2~T1 the log-normal dis- 
tribution does have a maximum. This implies that, depending on 
the pulsar flux, parallaxes can be biased towards larger as well as 
smaller values. 

Our software for estimation of the bias is available on-linfl 
both through a web interface and for download and off-line execu- 
tion. 



3 PARALLAX MEASUREMENT REVISIONS 

We have determined the bias-corrected confidence intervals for all 
known VLBI and timing parallax measurem ents found in the ATNF 
pulsar catalogue (Manches ter et alj|2005h . on Shami Chatterjee's 
dedicated webpag^fl and through NASA's Astrophysics Data Sys- 
tem Bibliographic Services. These results can be found in TablesQ] 
and[2] Note that in these tables and throughout the paper, we quote 

1 a error bars. 

In order to evaluate the potential of our bias-correction 
method, we compare historic measurements and their corrected val- 
ues to more precise, recent measurements in Section [3~Tl Next, we 
discuss the most significant revisions presented, in Section [3~/2l Pos- 
sible effects on models of the interstellar medium density and dis- 
tances derived from these models are discussed in Section [331 In 
Section [3~4l the potential impact of Lutz-Kelker bias on systems of 
particular interest for fundamental tests of gravitational theories is 
described. Finally, bias corrections for some non-radio pulsars are 
considered in Section [331 

3.1 Historic Parallax Measurements 

As discussed in Sj2. 1 l and shown by Equation[7] the size of the Lutz- 
Kelker bias is strongly dependent on the significance of the input 
parallax measurement. This implies that, as increasingly precise 
parallax measurements are published, those parallax values will 
converge towards decreasing bias. Comparing the historic paral- 
lax measurements collated in Table[2]to their most precise updated 

2 http://psrpop.phys.wvu.edu/LKbias 

3 Shami Chatterjee's pulsar parallax webpage can be found at: 
http://www.astro.cornell.edu/research/parallax 
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Figure 1. Example output from our Monte- Carlo simulations, ap plied to 
the PSR J1713+0747 parallax measurement of lCamilo et"afUl994l) (zoo = 
0.9 ± 0.3 mas). The top figure shows the volumetric and luminosity prior 
distributions (dashed and dot-dashed lines respectively) as well as their nor- 
malised product, the total prior distribution (solid line). The bottom figure 
shows the original measurement (dotted line), the prior distribution (dashed 
line) and their normalised product: the bias-corrected probability distribu- 
tion (solid line). Vertical lines show the peak and 1 a error bars of the mea- 
sured (dotted) and bias-corrected (solid line) probability distributions. 

values, we can now assess if these corrections have been towards 
the peak of the prior distribution, or away from it, respectively in- 
dicated by "T" and "A" in column 5 of Table Q] (The peak of the 
prior probability density function is determined from the pulsar ra- 
dio flux and the volumetric probability density function towards the 
pulsar and is provided SIS ZU~Prior in column 5 of TableQ]) 

This inspection shows that 31 out of 43 revisions were in the 
direction predicted by the prior distributions. The chance of this 
happening at random is less than 0.5% so this confirms the de- 
scribed bias effects are present in the data. It does not, however, 
imply bias-correction can resolve all these variations. Specifically, 
measurements with underestimated uncertainties will be under- 
corrected, as is supposedly the c ase for the PSR J0953+0755 m ea- 
surement bv lGwinn etaTI il9Sdt) (see §4 in lBrisken et"ai]|2000l . for 
a discussion of this measurement). We note that nearly half of 
all revisions listed in Table [2] are significant at the 1 a level, ir- 
respective of bias-correction. This strongly suggests more wide- 
spread underestimation of parallax uncertainties. Also, in cases 
where stronger systematic biases are present, application of the 
Lutz-Kelker correction will result in equally unreliable results. For 
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Table 1. Summary of published parallax values. This table contains the most up-to-date parallax measurements of pulsars derived from VLBI, pulsar timing 
and optical astrometry. Given are the pulsar name, published parallax value too. f ux at 1400MHz S140O: bias-cor rected parallax CTc or r , t he peak of the 
prior d istribution TOp r j or and any related references. Unless another reference is given, the fluxes were compiled from Lorimer et al. t 1995) and Kra mer et alj 
1 1998). Pulsars for which the parallax was determined optically are identified by asterisks next to their name. Any previous measurements are collated in Table 
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Pulsar 


name 




5l400 




rop r i or 


Ref. 


Jname 


Bname 


(mas) 


(mjy) 


(mas) 


(mas) 




J0030+0451 




3.3 ± 0.9 


0.6 


1.8+ ° 


0.5+2/6 

— U.2 


Lommen et al. (2006); Lommen et al. (2000) 


J0034-0721 


B0031-07 


93+°'°! 

yJ -^- J — 0.07 


1 1 


93+2 : 2? 

0.07 


1 o+J,- 3 


Chatterjee et al. (2009) 


J0108— 1431 




4.2 ± 1.4 


1.0 


i - u — 0.6 


9+2.1 


Deller et al. (2009) 


J0139+5814 


B0 136+57 


0.37 ± 0.04 


4.6 


0.37 ± 0.04 


38+2"?? 
UJO — 0.13 


Chatterjee et al. (2009) 


J0332+5434 


B0329+54 


0.94 ±0.11 


203 


0.93 ±0.11 


061 + 1 -1 
u - ul - 0.24 


Brisken et al. (2002) 


J0358+5413 


B0355+54 


0.91 ±0.16 


23 


0.87 ±0.16 


39+ ' 73 
u;)y -0.13 


Chatteriee et al. (2004) 


J0437-4715 




6.396 ± 0.054 


142 


0-»3_o .055 


1 01+ 1 ' 79 
lul -0.37 


Deller et al. (2008) 


J0452-1759 


B0450-18 


0.65+ ; 4 


5.3 


Q.69+°J? 


ftQ+0-88 
uoy -0.22 


Chatteriee et al. (2009) 


J0454+5543 


B0450+55 




13 




. in 7Q 

u#y -0.17 


Chatteriee et al. (2009) 


J0538+2817 






1.9 


o.»±S:m 


n O-7+0.39 
u - z/ -0.08 


Chatteriee et al. (2009); Lewandowski et al. (2004) 


T061 3— 0700 




80 + 35 


1.4 


046 +0.33 


31+0.40 

UJ1 -0.10 


Verbiest pt al f?00Q*l 


J0630— 2834 


B0628— 28 


3.0 ± 0.4 


23 


2.9 ± 0.4 


65+i'2i 

u.uj_ 22 


Deller et al. (2009) 


J0633+1746* 




6.4 ± 1.8 




17+2 -^f 

—0.01 


15+2'?? 
J — 0.01 


Caraveo et al. (1996) 


J0659+1414 


B0656+14 


3.47 ± 0.36 


3.7 


3 32+°'?S 

J - Jz '— 0.37 


"■^'-0.14 


Brisken et al. (2003) 


J0720— 3125* 




2.77 ± 0.89 




24+i'S 

"■ z -0.08 


20+2'?'? 
ui "-0.06 


Kaplan et al. (2007) 


J0737-3039A 




0.87 ±0.14 


1.6 




21+0- 29 
-0.06 


Deller et al. (2009); Bureav et al. (2006) 


J0751+1807 




1.6 ±0.8 


3.2 


6+ ' 7 
u.o_ 2 


045+0.59 
u/ "-0.14 


Nice et al. (2005) 


J08 14+7429 


B0809+74 


2.31 ±0.04 


10 


9 Tfl+0.05 
z - 5u -0.04 


77 +0.97 
u '"-0.26 


Brisken et al. (2002) 


J0820-1350 


B0818-13 


51+0 03 
u - ,1 -0.04 


7 


5 1 +°' 03 
u - ,1 -0.04 


40+0' 62 
u '^ u -0.13 


Chatteriee et al. (2009) 


J0826+2637 


B0823+26 


2.8 ± 0.6 


10 


2.5 ± 0.6 


7 c+0.99 
u ' /;, -0.24 


Gwinn et al. (1986) 


J0835— 4510* 


B0833— 45 


3.5 ± 0.2 


1100 


3.5 ± 0.2 


83+i'?i 


Dodson et al. (2003); Backer & Fisher (1974) 


J0922+0638 


B09 19+06 


0.83 ± 0.13 


4.2 


0.83 ± 0.13 


78+?,'*? 

w u — 0.25 


Chatteriee et al. (2001) 


T0953+0755 

JU7JJTU 1 


R0950+08 


182 + n 07 


84 


3 81 + 07 


1 20+ 1 ' 50 
lzu -0.42 


Rrisken et al OWO\ 

UllAiYCll til. ^ i-UVli J 


J1012+5307 




1.22 ± 0.26 


3 


1 09+2 -11 

'■"'-0.23 


65+2'o? 

o.uj_o 21 


Lazaridis et al. (2009) 


ti 022+1 001 




1.8 ± 0.3 


3 


1.7 ± 0.3 


62+ ' 74 

uoz -0.19 


Verhiesf pt al l'?00Q'l 

vui uirjjL til. \sL\j\jyj 


J1024-0719 




1.9 ±0.8 


0.66 


07 +0.7 


045+0.50 
u - #:) -0.13 


Hotan et al. (2006) 


J 1045 -4509 




3.3 ± 1.9 


3 


025 +0.47 
u - z3 -0.08 


u ' z ^-0.07 


Verbiest et al. (2009) 


J1136+1551 


B1133+16 


2.80 ±0.16 


32 


n Q+0.15 


1 90+1. 40 
lzy -0.44 


Brisken et al. (2002) 


J1239+2453 


B 1237+25 


1.16 ± 0.08 


10 


1 1 -T+0 07 
la '-0.09 


, , _+l 20 

!- 17 Iai9 


Brisken et al. (2002) 


J 1300+ 1240 


B1257+12 


1.3 ±0.4 


2 


, i_0 A 

1 0^ 


_ „ , _i_n 79 
0.64i°;™ 


Wolszczan et al. (2000) 


J1456-6843 


B1451-68 


2.2 ± 0.3 


80 


2.1 ± 0.3 


37+2'?2 
^•^ ' —0.13 


Bailes et al. (1990); Manchester et al. (1980) 


J1509+5531 


B 1508+55 


0.47 ± 0.03 


8 


0.48 ± 0.03 


82+i'2? 

"■o^_0.25 


Chatterjee et al. (2009) 


J1537+1 155 


B1534+12 


0.98 ± 0.05 


0.6 


0.97 ± 0.05 


O60+2' 59 , 

u.uo_ lg 


Stairs et al C20021 

*JLtlllCJ w L til. \ £*\S\J£, f 


J 1543+0929 


B 1541+09 


0.13 ± 0.02 


5.9 


0.16 ± 0.02 


75+2'2? 

u - ' — 0.24 


Chatteriee et al. (2009") 


T1SSQ AdW 

j ujy +^t_)o 


pU55fi A A 




40 


3(SK+ 081 
U.JO»_ 07g 


U - ZZ -0.07 


DpIW pt nl fOOTIQV Frmmlrmt pt nl MQQ?^ 
UCllCl cL til. (ZUUyj, rUHlalUilL cL til. \ lyyZ.) 


J1600-3053 




0.20 ±0.15 


3.2 


02 ,+0.12 
u - zl -0.06 


2] +0.32 
u - z, -0.06 


Verbiest et al. (2009); Jacobv et al. (2007) 


J1643-1224 




2.2 ± 0.4 


4.8 


1.9 ± 0.4 


90+ ' 34 
u - zu -0.05 


Verbiest et al. (2009) 


J1713+0747 




0.94 ± 0.05 


8 


0.93 ± 0.05 


2 ,+0.42 
UZJ -0.07 


Verbiest et al. (2009) 


J1744— 1134 




2.4 ±0.1 


3 


2.4 ±0.1 


14+0 18 
ul4 -0.03 


Verbiest et al. (2009) 


J1856-3754* 




6.2 ± 0.6 




6.0 ± 0.6 


17+0- 20 
u - 1/ -0.04 


van Kerkwiik & Kaplan (2007) 


J 1857+0943 


B 1855+09 


1.1 ±0.2 


5 


0.9 ± 0.2 


13+0 16 


Verbiest et al. (2009) 


J1900-2600 


B1857-26 


0.5 ± 0.6 


13 


n ,o+0.30 
U - 1S -0.05 


17+0- 35 
ua/ -0.04 


Fomalont et al. (1999) 


J1909-3744 




0.79 ± 0.02 


3 


0.79 ± 0.02 


1Q+0.30 
Uly -0.05 


Verbiest et al. (2009); Jacobv et al. (2003) 


J1932+1059 


B1929+10 


2.77 ± 0.07 


36 


2.76 ± 0.07 


19+ ' 44 


Chatterjee et al. (2004); Hobbs et al. (2004) 


J1935+1616 


B1933+16 


u - zz -0.12 


42 


70+0- 33 
u - zu -0.06 


1Q+0- 46 


Chatteriee et al. (2009) 


J1939+2134 


B 1937+21 


0.13 ±0.07 


10 


14+005 
u - 14 -0.03 


,5+0.25 
u - 1:> -0.04 


Verbiest et al. (2009) 


J20 18+2839 


B2016+28 


1.0 ±0.1 


30 


1.0 ±0.1 


22+0- 49 
u - zz -0.07 


Brisken et al. (2002) 


J2022+2854 


B2020+28 


0.37 ±0.12 


38 


O34+O.12 
U - J *-0.10 


fl na+0.53 
u - ZJ -0.07 


Brisken et al. (2002); Gomez-Gonzalez & Guelin (1974) 


J2022+5154 


B2021+51 


0.50 ± 0.07 


27 


0.50 ± 0.07 


39+ ' 74 
UJy -0.13 


Brisken et al. (2002) 


J2048-1616 


B2045-16 


1 05+ 03 
i- u J_0.02 


13 


1 05+ 03 


ifi+0- 85 
u; > D -0.18 


Chatteriee et al. (2009) 
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Table 1 - continued 



Pulsar name vjq Si400 ^Corr ro Prior Ref- 

Jname Bname (mas) (mJy) (mas) (mas) 



J2055+3630 


B2053+36 


0.17 ±0.03 


2.6 


0.18 ± 0.03 


77+0.39 
u - z/ -0.08 


Chatteriee et al. (2009) 


J2124-3358 




3.1 ±0.6 


1.6 


2.7 ± 0.6 


0A2±°f 3 


Verbiest et al. (2009) 


J2129-5721 




1.9 ±0.9 


1.4 


4+0.7 

uq -0.2 


0.36lg;« 


Verbiest et al. (2009) 


J2144-3933 




6.05 ± 0.56 


0.8 






Deller et al. (2009) 


J2145-0750 




1.6 ±0.3 


8 


1.5 ± 0.3 


52+ ' 72 


Verbiest et al. (2009) 


J2157+4017 


B2154+40 


0.28 ± 0.06 


17 


uzy -0.05 


4 +0.70 


Chatteriee et al. (2009) 


J23 1344253 


B2310+42 


uy ->-0.07 


15 


uyz -0.07 


56 +0.85 

u.->o_ 19 


Chatteriee et al. (2009) 



PSR J0437-4715 for example. Iverbiest et al] d2008h have demon- 
strated that the parallax value published by van Straten et al] ( 200 1 ) 
differed from later measurements because of inaccuracies in the So- 
lar System Ephemeris models used. 



3.2 Significant Parallax Revisions 

As follows from Table Q] out of the 56 pulsars that currently have 
parallax measurements from radio or optical data, only 9 have 
a bias larger than one standard deviation. In the case of PSRs 
J0633+1746 and J0720-3125 this is mostly due to the absence of 
intrinsic luminosity information, as described in Section [3~5l while 
for PSR J075 1+1807 the parallax bias is mostly caused by the low 
significance of the measurement. For the remaining 6 pulsars the 
significance of the bias is caused by the fact that the luminosity 
and volumetric prior peak on the same side of the measured value. 
This causes the two components to the prior to reinforce each other, 
resulting in more significant bias corrections. 

The bias correction of PSR J0108— 143 1 is of particular inter- 
est, since this pulsar was on ce assumed to be the closest known neu- 
tron star jTauris et all 19941) . based on its DM distance of D1994 = 
130 pc as derived fro m the Galactic electron d ensity model of 
ITavlor & Corded dl993h . The updated model of ICordes & Laziol 
d2002h i ncreased the distanc e to D2002 = 200 pc and a recent VLBI 
parallax dDeller e t al. 2009) further increased this value to D2009 = 
240l"e;i 4 pc. Our bias-corrected parallax rocorr = l-6t"J'g places 



625 tl 7 A pc, strong! 



the pulsar at a distance of 7Jjoio8-i43 

suggesting that the 3 a parallax measurement o flDeller et al] d200' 
is still underestimating the distance. 

The bias-corrected parallax of PSR B1541+09 is 1.5 a larger 
than the VLBI value. This pulsar has a Galactic latitude of 46°, 
which in co mbination with the par allax measurement of zuq = 
0.13 ± 0.02 dChatteriee et al]|2009l) would place it at a height of 
5.5 kpc above the Galactic disk. This is a highly unlikely height 
given the scale height of 330 pc for normal pulsars, suggesting that 
either the parallax measurement is strongly underestimated because 
of bias effects or systematic errors, or the scale height in our Galac- 
tic model is underestimated near the Galactic bulge. 



3.3 ISM Density Models 

As mentioned in Section [T] models of the Galacti c electron distri- 
butio n are based on pulsar parallax measurements ( Cordes & Laziol 
2002). Since the estimated biases for most measurements are small 
(as follows from Table[TJ and the uncertainties in the electron den- 
sity models are relatively large, bias correction is not expected to 



significantly affect such models. In future modelling efforts, how- 
ever, bias-corrected parallax values should result in more reliable 
models. 

Conversely, Galactic electron density models are used to es- 
timate pulsar distances based on the DM. While these distance es- 
timates are not trigonometric, the same reasoning as presented in 
Section [2~71 holds and DM-derived distances can be expected to be 
affected by Lutz-Kelker bias. In order to investigate the presence 
of Lutz-Kelker bias in DM distances, we consider all pulsars with 
trigonometric parallaxes that were not contained in the most recent 
Galactic electron density model and that have m easurement uncer- 
tainties lower than 20% dCordes & Lazioll2003l. their Figure 12). 
The DM distances and parallax values for these 21 pulsars are com- 
piled in Table[3] assuming a 20% uncertainty on the DM distance. 
Comparison of the DM distances with the more precise parallax 
measurements from pulsar timing and VLBI compiled in Table Q] 
shows that in 15 of the 21 cases the more precise measurement 
lies closer to the peak of the prior distribution cap r ior, indicating 
the presence of Lutz-Kelker bias. Bias-corrected DM-derived par- 
allaxes are also provided, though interpretation of these values is 
complicated by the fact that the uncertainties of DM distances are 
rough estimates. This analysis demonstrates that DM distances are 
likely to be underestimated, a fact that will need to be taken into 
account in research that makes use of such distances, for examp le 
in the determination of pulsar velocities as in iHobbs etal]d2005l) . 



3.4 Binary Pulsars 

One particular application where bias correction may prove cru- 
cial in future research, is the determination of gravitational wave 
emission from binary pulsar systems. This emission is measured 
through orbital period dec ay, predicted by general relativity to be 
dTavlor & Weisberdll982h : 



1927T f2nT Q 
~5~ V Ph 



5/3 



M c Mr> 1 + ie 2 + §e 4 



AL 



1/3 



(1 



\7/2 



with Pb the first derivative of the orbital period, Pb, Tq = 
4.925490947/iS, M p the pulsar mass, M c the companion mass, 
Mtot = M p + M c the total system mass and e the orbital ec- 
centricity. Given an independent measurement of M c and M p , a 
measurement of Pb can be u sed to test gravitational wave emis- 
sion theories, as first done by ITavlor & Weisber ^ dl982l) . Alterna- 
tivel y, it can be used to determine the pulsar and companion masses, 
as in iNice et al] d20ol> . Depending on the pulsar's proper motion 
and position in the Galaxy, however, certain kinematic terms may 
contaminate the Pb mea surement, as first described in detail by 
iDamour & Tavlorl dl99ll) . They showed that three extra terms need 
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Table 2. Historic parallax measurements of pulsars. Given are the pulsar name, the published parallax value vjq, the significance of the measurement roo/cr, 
the change of the measurement with respect to the previously published value Aroo, T/A indicating whether the measurement was corrected Towards or Away 
from rop r i or (listed in Table [TJ when compared to the most precise measurement available, the bias-corrected parallax value rocorr, its difference with the 
previous bias-corrected value Arocorr and the relevant publication. 



Pulsar 


TOO 


Significance 


Aroo 


Change w.r.t. 


TOQorr 


AaJCorr 


Ref. 


name 


(mas) 




(mas) 


ro Prior 


(mas) 


(mas) 




J0332+5434 


1.3 ± 1.7 


0.8 


- 


T 


7+ 1 - 1 


- 


Salter et al. (1979) 


(B0329+54) 


0.91 ±0.11 


8.3 


-0.39 ± 1.70 


- 


0.91 ±0.11 


0.21 ± 1.11 


Brisken et al. (2002) 


J0437-4715 


5.6 ±0.8 


7.0 


- 


A 


5.3 ±0.8 


— 


Sandhu et al. (1997) 




7.19 ± 0.14 


51.4 


1.59 ± 0.81 


T 


7.16 ± 0.14 


1.86 ± 0.81 


van Straten et al. (2001) 




o.j ±0.1 


OJ.O 


n on _i_ a 1 7 

— u.oy ± U.l f 


A 


6. J ±0.1 


— U.ob ± U.l ( 


riotan et al. (ZUUo) 




o.oj ± 0.31 


13.0 


r\ o c i f \ c n 

U.jj ± U.;>2 


T 


<fi +0.48 

6.395±g;g|| 


U.25 ± O.jj 


Verbiest et al. (2UUo) 




6.396± 0.054 


118.4 


-0.254 ± 0.513 


— 


-0.185 ± 0.543 


Deller et al. (2008) 


J0538+2817 


0.68 ±0.15 


4.5 


- 


A 




- 


Ne et al. (2007) 




u - /z -0.09 


6.0 


0.04 ±0.17 


- 


0-68iao9 


0.08 ± 0.17 


Chatteriee et al. (2009) 


J0613-0200 


2.1 ±0.6 


3.5 




T 


i-i±o:l 




Hotan et al. (2006) 




U.OVJ HI U.JJ 


1 3 

Z.J 


i q _l n 7 




46+ 03 


n (\a -i- n 

— U.\J*± HZ u.uo 


veiDiesi ei di. \zx)\}y) 


J0737-3039A 


3±2 


1.5 




T 


0.21+gg5 




Kramer et al. (2006) 




0.87 ±0.14 


6.2 


-2.13 ± 2.00 


_ 


0.80+°^ 


0.59 ± 0.38 


Deller et al. (2009) 


J0814+7429 


1.8 ± 7.1 


0.3 


_ 


A 




_ 


Salter et al. (1979) 


(B0809+74) 


2.31 ±0.04 


57.8 


0.51 ±7.10 


_ 


2.31 ±0.04 


1.61 ± 1.00 


Brisken et al. (2002) 


J0835-4510 


3.4 ± 0.7 


4.9 


_ 


A 


3.1 ±0.7 


_ 


Caraveo et al. (20011 


(B0833-45) 


3.5 ± 0.2 


17.5 


0.1 ±0.7 


- 


3.5 ± 0.2 


0.4 ± 0.7 


Dodson et al. (2003) 


J0922+0638 


0.31 ±0.14 


2.2 


_ 


T 


0.43l°.f 




Fomalont et al. (19991 


(B09 19+06) 


0.83 ± 0.13 


6.4 


0.52 ±0.19 




o.83+° ; " 

— U. lo 


0.4 ± 0.2 


Chatteriee et al. (20011 


J0953+0755 


8.3 ± 8.8 


0.9 




T 


L '^— 0.4 




Salter et al. (1979) 


(B0950+08) 


7.9 ±0.8 


9.9 


-0.4 ±8.8 


T 


7.8+^ 
— o.y 


6.6 ± 2.3 


Gwinn et al. (1986) 




3.6 ± 0.3 


12.0 


—4.3 ± 0.9 


A 


3.5 ± 0.3 


—4.3 ± 1.0 


Brisken et al. (2000) 




3.82 ± 0.07 


54.6 


0.22 ± 0.31 




3.81 ± 0.07 


0.31 ± 0.31 


Brisken et al. (2002) 


J1022+1001 


3.3 ± 0.8 


4.1 




T 


2 4+°-^ 

^■^-0.8 




Hotan et al. (2004) 




Z.J zn u.t 


O.J 


n a 4- n q 
— u.o m u. y 


T 


24 +0.4 


4-1 03 

U.U ZL1 1 .U J 


nuidii ei ai. 




1 .0 ZH U.J 


6 


7 4-0 ^ 
— u. i m u.o 




j 7 z_ o 3 


7 4- ^8 


veiuiesi ci di. \jLU\J7) 


J1239+2453 


-6.2±5.2 


-1.2 


- 


T 


1 o +1 -° 

i * u — 3 


- 


Salter et al. (1979) 


(B 1237+25) 


1.16 ±0.08 


14.5 


7.36 ± 5.20 




j 17 +o:o7 


0.17 ± 1.00 


Brisken et al. (2002) 


J 1509+5531 


0.415± 0.037 


11.2 




T 


0.430+S 




Chatteriee et al. (2005) 


(B1508+55) 


0.47 ± 0.03 


15.7 


0.055 ± 0.048 


_ 


0.48 ± 0.03 


0.050 ± 0.045 


Chatteriee et al. (2009) 


J1537+1155 


0.91+2?? 

— 0.14 


6.5 




A 


0.89+0?? 

" " v — 0.13 




Stairs et al. (1998) 


(B 1534+12) 


0.98 ± 0.05 


19.6 


0.07 ± 0.05 




0.97 ± 0.05 


0.08 ± 0.05 


Stairs et al. (2002) 


J1713+0747 


0.9 ± 0.3 


3.0 




A 


6+S-^ 

"■"-9.2 




Camilo et al. (1994) 




89 + 08 

W.t3/ _1_ U.UO 


11.1 


— 01 + 31 


A 


87 + 08 


0.27 ± 0.31 


Splaver et al. (2005) 




1 10 ± 05 


22.0 


0.21 ± 0.09 


T 


1 1 ± 05 


0.23 ± 0.09 


Hotan et al. (2006) 




0.95+°-°t 
— u.uo 


15.8 


-0.15 ± 0.08 


T 


0.94+ -08 

— O.OD 


-0.16 ±0.08 


Chatteriee et al. (2009) 




0.94 ± 0.05 


18.8 


—0.01 ± 0.07 




0.93 ± 0.05 


—0.01 ± 0.07 


Verbiest et al. (2009) 


J1744-1134 


2.8 ± 0.3 


9.3 




T 


2.6 ± 0.3 




Toscanoetal. (1999) 




2.1 ± 0.2 


10.5 


-0.7 ±0.4 


A 


2.0 ± 0.2 


-0.6 ±0.4 


Hotan et al. (2006) 




2.4 ±0.1 


24.0 


0.3 ±0.2 




2.4 ±0.1 


0.4 ± 0.2 


Verbiest et al. (2009) 


J1856-3754 


16.5 ± 2.3 


7.2 




T 


15.2 ±2.4 




Walter (2001) 




7 ± 2 


3.5 


-9.5 ± 3.1 


T 


17+ - 57 

vl 1 -0.04 


-15.0 ±2.5 


Kaplan et al. (2002) 




8.5 ± 0.9 


9.4 


1.5 ±2.2 


T 


7.2 ± 1.6 


7.03 ± 1.70 


Walter & Lattimer (2002) 




6.2 ± 0.6 


10.3 


-0.8 ± 1.1 




6.0 ± 0.6 


5.83 ± 1.71 


van Kerkwiik & Kaplan (2007) 


Jl 857+0943 


1.2 ±0.5 


2.4 




T 


n 14+036 
ua4 -0.04 




Rvba & Taylor (1991) 


(B 1855+09) 


1.1 ±0.3 


3.7 


-0.1 ± 0.6 


A 


0.7 ± 0.4 


0.56 ± 0.54 


Kasni et al. (1994) 




1.1 ±0.2 


5.5 


0.0 ± 0.4 




0.9 ± 0.2 


0.2 ± 0.5 


Verbiest et al. (2009) 


J1909-3744 


0.88 ± 0.03 


29.3 




T 


0.88 ± 0.03 




Jacobv et al. (2005) 




0.88 ± 0.02 


44.0 


0.0 ± 0.04 


T 


0.88 ± 0.02 


0.0 ± 0.04 


Hotan et al. (2006) 




0.79 ± 0.02 


39.5 


-0.09 ± 0.03 




0.79 ± 0.02 


-0.09 ±0.03 


Verbiest et al. (2009) 


J1932+1059 


21. 5± 8.0 


2.7 




T 


n iq+0.58 
U - ly -0.06 




Salter et al. (1979) 


(B1929+10) 


5.0± 1.5 


3.3 


-16.5 ±8.1 


T 


n 99+2.8 


0.03 ± 0.58 


Camilo (1995) 




3.02± 0.09 


33.6 


-1.98 ± 1.50 


T 


3.01 ± 0.09 


2.79 ± 2.80 


Brisken et al. (2002) 




2.77± 0.07 


39.6 


-0.25 ± 0.11 




2.76 ± 0.07 


-0.25 ±0.11 


Chatteriee et al. (2004) 
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Table 2 - continued 



Pulsar 


— 1 1 


Significance 


A777Q 


Change w.r.t. 


name 


(mas) 


vjq/o- 


(mas) 


ro Prior 


J1939+2134 


0.12 ± 0.08 


1.5 


- 


T 


(B 1937+21) 


0.13 ± 0.065 


2.0 


0.01 ±0.10 


- 


J20 18+2839 


0.9± 3.6 


0.3 


- 


A 


(B2016+28) 


1.03± 0.10 


10.3 


0.13 ± 3.60 




J2022+2854 


-1.3±4.6 


-0.3 




T 


(B2020+28) 


0.37± 0.12 


3.1 


1.67 ±4.60 




J2022+5154 


1.8± 4.9 


0.4 




T 


(B2021+51) 


0.95± 0.37 


2.6 


-0.85 ±4.91 


T 




0.50± 0.07 


7.1 


-0.45 ±0.38 




J2048-1616 


1.71± 0.91 


1.9 




T 


(B2045-16) 


± - uo -0.02 


35.0 


-0.66 ±0.91 




J2124-3358 


4±2 


2.0 




T 




3.1 ± 0.55 


5.6 


-0.9± 2.1 




J2145-0750 


2.0 ±0.6 


3.3 




T 




1.6 ±0.25 


6.4 


-0.4 ±0.7 





(mas) 



(mas) 



Ref. 



H1.06 
-0.04 
,,+0.05 
-0.03 
■,+0.48 
-0.07 



0.14Z 
0.14;! 
0.22^ 
1.00 ±0. id 
0.23j 
0.34^ 
0.4j 

o.7r 



5+0.50 
-0.07 
,,+0.12 
-0.10 
,+0.8 
-0.1 
(-0.37 
-0.26 
0.50 ± 0.07 
5+0.82 
-0.26 
,+0.03 
-0.02 
^+0.8 
-0.2 
7+0.56 
-0.61 
!+0.6 
-0.5 

1.5 ±0.3 



0.73: 
1.05"* 

0.5j 
2.67"' 

1.4 4 



0.0 ± 0.08 

0.78 ± 0.49 

0.11 ±0.51 

0.31 ±0.84 
-0.21 ±0.27 

0.32 ± 0.82 

2.17 ± 1.01 

0.1 ±0.7 



Kaspi et al. (1994) 
Verbiest et al. (2009) 

Salter et al. (1979) 
Brisk en et al. (2002) 

Salter et al. (1979) 
Brisk en et al. (20021 

Salter et al. (1979) 
Campb ell et al. (1996) 
Brisken et al. (2002) 

Deller et al. (2009) 
Chatteri ee et al. (2009) 

Hotan et al. (2006) 
Verbiest et al. (2009) 
Lohmer et al. (2004) 
Verbie st et al. (2009) 



Table 3. Lutz-Kelker corrections for DM- derived distan ce estimates. Given are the pulsar name, the DM-derived distance according to the NE2001 model for 
Galactic electron density Ddm ICordes & Lazidl200ll) . the parallax edm derived from D, assuming a 20% uncertainty in the DM-distance and the Lutz- 
Kelker corrected parallax ■cuum Corr- This table only contains distances for pulsars that have accurate (relative errors less than 20%) independent distances 
(from timing or VLBI), which were obtained after the NE2001 model was created. 



Pulsar 


£>DM 




^DM.Corr 


name 


(kpc) 


(mas) 


(mas) 


J0034-0721 


0.41 


2.4(5) 


2.2 ± 0.5 


J0139+5814 


2.88 


0.35(7) 


35+ ' 07 
u -"-0.06 


J0437-4715 


0.14 


7(1) 


6 3+ 14 
°- 3 —1.6 


J0454+5543 


0.67 


1.5(3) 


1.4 ±0.3 


J0538+2817 


1.22 


0.8(2) 


074 +o.i6 


J0737-3039A 


0.52 


1.9(4) 


1.6 ±0.4 


J0820-1350 


1.99 


0.5(1) 


u -^ y -0.09 


J0835-4510 


0.24 


4.2(8) 


, n+0.8 
J - y -0.9 


J1012+5307 


0.41 


2.4(5) 


2.1 ± 0.5 


J1022+1001 


0.45 


2.2(4) 


1.89 ± 0.45 


J1509+5531 


0.99 


1.0(2) 


1.0 ±0.2 


J1543+0929 


3.49 


0.29(6) 


0.33 ±0.05 


J17 13+0747 


0.89 


1.1(2) 


1.0 ±0.2 


J1744-1134 


0.41 


2.4(5) 


1 q+0.5 


J1909-3744 


0.46 


2.2(4) 


1.8 ±0.5 


J1932+1059 


0.34 


2.9(6) 


2.5 ± 0.6 


J2048-1616 


0.56 


1.8(4) 


1.6 ±0.4 


J2055+3630 


4.62 


0.22(4) 


0.22 ± 0.04 


J2124-3358 


0.27 


3.7(7) 


i o+0.7 
J - z -0.9 


J2145-0750 


0.57 


1.8(4) 


1.6 ±0.4 


J23 13+4253 


1.25 


0.8(2) 


0.78 ± 0.15 



to be considered and corrected for: t he Shklovskii ter m based on 
the pulsar proper motion and distance ( Shklovskii 1970); the appar- 
ent acceleration caused by the gravitational potential of the Galaxy, 
dependent on the pulsar's Galactic latitude and height above the 
Galactic plane; and the apparent acceleration caused by differen- 
tial Galactic rotation, which is also dependent on the Galactic lat- 
itude and therefore the distance of the pulsar. Since all three non- 
relativistic contributions to the measured orbital period derivative 
Pb are dependent on the pulsar distance, accurate estimation of 
these contaminating factors and subsequent determination of grav- 



itational wave emission or pulsar mass estimates, is dependent on 
an accurate parallax determination. 

Presently, the orbital decay Pb has only been determined for 
16 binary pulsar systems, of which only four have a parallax mea- 
surement. In the cases of PSRs J0437-4715 and J1909-3744 
the parallax measurement is precise enough to remain unaffected 
by bias effects and in the cases of PSRs J0737-3039A/B and 
J0751+1807 the uncertainty in the measured Pb value is substan- 
tially larger than the expected contributions from kinematic and 
Galactic effect dKramer et alj2006l;lNic"e et alj2005t) . 
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3.5 Optical Pulsar Parallax Measurements 

A total of four neutron stars have had their parallaxes measured 
though optical observations with the Hubble Space Telescope 
(PSRs J0633+1746, J0720-3125, J0835-4510 and J1856-3754). 
Of these four, only the Vela pulsar (PSR J0835-4510) has been 
detected at radio wavelengths. Since an optical intrinsic luminos- 
ity function of neutron stars has thus far not been determined, we 
disregard the luminosity prior for the remaining three optical neu- 
tron stars. Because surveys for radio pulsars are generally biased 
towards bright pulsars, most known radio pulsars lie at the bright 
end of the intrinsic luminosity function. This implies that the lu- 
minosity prior tends to counterbalance the volumetric prior, which 
is biased towards small parallax values (as seen in Figure [T}. The 
exclusion of the luminosity prior for non-radio neutron stars im- 
plies that the effect of the volumetric prior becomes much stronger, 
which means higher measurement precision is required before any- 
thing definite can be said about the bias in the neutron star distance. 
A clear exampl e of this is PS R J1856— 3754 for which the 3.5 a 
measurement of Kapl an et alj J2002h is insufficient to overcome the 
steep volumetric prior. For pulsars with luminosity information, a 
3.5 a measurement generally does suffice to restrict the bias within 
the 68% confide nce interval as illus trated by the PSR J0613-0200 
measurement of Hot an et al] d2006l) . We note that the most precise 
parallaxes for PSRs J0633+1746 and J0720-3125 are only 3.6 a 
and 3.1 a measurements respectively, which is insufficient to claim 
confidence in the inferred distance because of the poor estimate of 
the bias. 

X-ray luminosities of neutron stars are commonly used to de- 
rive neutron star radii which in turn are used to place constraints 
on equations of state for dense nuclear matter. One of the limiting 
factors in such analyses are the high distance uncertainties since 
distances are mainly derived from mostly inprecise methods such 
as interstellar dispersion, supernova remnant associations and inter- 
stellar ab sorption lines, wh ile only few parallax measurements are 
available IPage et al] {2004). This lack of precise distances causes 
even 2 and 3 a parallax measurements to carry significant weight 
in these analyses. As our work shows, without prior information on 
the intrinsic luminosity distribution of these sources, it cannot be 
confidently claimed that parallax measurements with such low sig- 
nificance are accurate enough to be used because they are likely to 
underestimate the distance to the neutron star. Based on equ ations 
56 and 58 as well as Figure 2A o flLattimer & Prakasbl d2007). these 
underestimated distances can be expected to result in underestima- 
tions of neutron star masses and radii, slightly biasing the analysis 
to prefer softer equations of state. 



4 CONCLUSIONS 

We ha ve reanalysed bias effects first discussed by iLutz &~K elker 
dl973h and described Monte-Carlo simulations that aim to cor- 
rect for these biases. Comparison of historic parallax values to the 
most recent measurements confirms that the bias effects are present 
in observations, though the correction is complicated by system- 
atic measurement errors and underestimation of parallax uncertain- 
ties. Correction for the described biases may improve some pulsar- 
timing tests of gravitational wave emission and may slightly influ- 
ence future density models of the interstellar medium. Finally, we 
conclude that optical parallax measurements should be used with 
caution since we cannot acurately quantify the possible bias in the 
inferred distance estimates. 
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